Oscillatory Curie temperature of two-dimensional ferromagnets 
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The effective exchange interactions of magnetic overlayers Fe/Cu(001) and Co/Cu(001) covered by 
a Cu-cap layer of varying thickness were calculated in real space from first principles. The effective 
two-dimensional Heisenberg Hamiltonian was constructed and used to estimate magnon dispersion 
laws, spin-wave stiffness constants, and overlayer Curie temperatures within the mean-field and 
random-phase approximations. Overlayer Curie temperature oscillates as a function of the cap- 
layer thickness in a qualitative agreement with a recent experiment. 
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The Curie temperature is one of the most important 
characteristics of ferromagnets. In particular, the Curie 
temperature of low-dimensional systems such as ultrathin 
films is of considerable interest. In a recent experimental 
study, Vollmer et al. |ij have shown that (i) the Curie 
temperature of fcc(001)-Fe ultrathin films on a Cu(001) 
substrate is considerably modified upon coverage by a 
Cu-cap layer, and (ii) that it varies in a non-monotonous 
manner as a function of the Cu cap layer thickness, which 
indicates an oscillatory variation. An oscillatory behav- 
ior of the Curie temperature as a function of the spacer 
thickness was also found for fcc(001)-Co/Cu/Ni trilay- 
ers Q|. Such a behavior clearly cannot be explained 
within a localized picture of magnetism and calls for a 
first-principles theory of the Curie temperature in itiner- 
ant ferromagnets. In spite of considerable efforts in last 
decades a first-principles calculation of the Curie temper- 
ature in the framework of itinerant magnetism, in partic- 
ular for low-dimensional systems, remains a very serious 
challenge. 

One therefore has to rely upon some approximation 
schemes in order to calculate the Curie temperature of 
itinerant ferromagnets. A particularly simple and yet 
accurate approach consists in a mapping of the compli- 
cated itinerant electron system onto an effective Heisen- 
berg model (EHM), H = —J^i^jJij e i ' e ji where 
and ej are the unit vectors of the magnetic moments 
at sites i and j, and the effective exchange interactions 
(EEIs) Jij between any pair of magnetic moments are de- 
termined from first-principles |^-^| ■ The thermodynamic 
properties of the ferromagnet including determination of 
the Curie temperature can be then calculated from the 
EHM by using statistical mechanical methods. A simple 
mean field approximation (MFA) fails in many cases due 



to its neglect of collective excitations (spin- waves) , and 
more sophisticated approximations, such as the Green 
function method within the random phase approxima- 
tion (GF-RPA) p0[ , are preferable. The success of this 
two-step approach relies upon the fact that it provides 
an almost exact description of low-lying magnetic exci- 
tations (spin-waves) which give the largest contribution 
to the Curie temperature. On the other hand this ap- 
proach completely disregards longitudinal fluctuations of 
magnetic moments such as the Stoner excitations and it 
therefore is not suitable to describe ferromagnets with 
small exchange splitting such as, e.g., fcc-Ni, in which 
exist Stoner excitations with a rather low energy. We 
have recently applied this approach to bulk bcc-Fe, fcc- 
Co, and fcc-Ni and obtained a reasonable agreement with 
experimental Curie temperatures of Fe and Co, but not 
for Ni H similarly as in recent calculations based on 
the adiabatic spin-wave theory [Q or an alternative first- 
principles theory of spin fluctuations based on idea of 
generalized Onsager cavity field |pd| . 

In the present paper we wish to calculate exchange in- 
teraction parameters, spin-wave stiffness constants, and 
Curie temperatures of two-dimensional monolayers of Fe 
and Co and, in particular, to investigate the influence of 
the substrate and of the cap layer. We find that (i) the 
exchange parameters, spin-wave stiffness constants, and 
Curie temperatures are strongly modified by the presence 
of a metallic substrate and/or a cap layer, and that (ii) 
they exhibit an oscillatory variation with the cap layer 
thickness. This behavior is due to the Ruderman-Kittel- 
Kasuya-Yoshida (RKKY) character of exchange interac- 
tions in itinerant ferromagnets. Our results are in a good 
qualitative agreement with the observations of Vollmer 
et al. |IJ for which they provide the most natural expla- 
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nation. The same theory can be used to interprete the 
experiment of Ney et al. 0], but the detailed analysis 
deserves a separate study. 

The electronic structure of the system was determined 
in the framework of the first principles tight-binding lin- 
ear muffin-tin orbital method (TB-LMTO) generalized 
to surfaces JT3] . In the framework of the magnetic force- 
theorem [| 13 1, the expression for the EEIs between two 
sites i and j anywhere in the system is p|fL4[ 



i- / Imtr L {^(z)4.(z)^(z)4(z)} dz . (1) 



Jij — 



Here tr^ denotes the trace over the angular momentum 
L = (£m), 5i(z) = P}(z) - Pj-(z) where Ff(z) are L- 
diagonal matrices of potential functions of the TB-LMTO 
method (a =f, j), energy integration is performed in the 
upper half of the complex energy plane over a contour 
C starting below the bottom of the valence band and 
ending at the Fermi energy, and gfj(z) are the site off- 
diagonal blocks of the system Green function correspond- 
ing to a given geometry. Possible lattice and/or layer 
relaxations at the overlayer are neglected. The inter- 
site Green functions gfj(z) can be evaluated either in the 
real space by using the cluster approach |3| , the recursion 
method ||,[l5||, or, as it is done in the present paper and in 
Ref. , by the Bloch transformation which employs the 
two-dimensional translational symmetry of a given layer 
(for more details concerning the computational method 
see Ref. JIJ] ) . We have calculated the EEI pairs J,j up to 
101-shells of the fcc(OOl) surface (i.e., up to the distance 
of about 10 a, where a is the lattice constant of the fee 
lattice). Such a large number of the EEIs is needed, in 
particular, for an accurate estimate of the spin-wave stiff- 
ness constant in a real space, as it is also known for the 
bulk case 0,^1 ■ In actual calculations the sites i,j were 
limited to the magnetic layer, which is a good approxi- 
mation in view of the smallness of the moments induced 
in the Cu. The spin-wave spectrum -E(q|| ), the spin- wave 
stiffness constant D, and the Curie temperatures T c MFj4 
and T PPA are expressed, respectively, in terms of the 
EEIs as follows 

£(qn) = ^ J °* ( x - cx pNi ' R *)) + A > 
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Here, Nu is the number of sites per layer, the qy-sum 
extends over the fcc(OOl) surface Brillouin zone, [Ib is 
the Bohr magneton, R oi = |R — Rj is the interatomic 
distance, M is the magnetic moment per atom, and A 
is the magnetic anisotropy energy. It should be noted 



TABLE I. Calculated values of effective exchange interac- 
tions for the first 10 shells of fcc(OOl) Fe- and Co-magnetic 
layers: a free-standing (fs) layer, an overlayer (ov) on 
fcc(001)-Cu, and an embedded (em) layer in fcc-Cu. Num- 
bers of atoms in a given shell and the corresponding shell 
radii (in units of lattice constants) are given in parenthesis 
and square brackets, respectively. Corresponding values of 
magnetic moments M are also given. 
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that T^ FA can be evaluated directly by using the one- 
site rotation term Jo (expressed in terms of the site- 
diagonal element of the magnetic layer Green function 
similarly as its bulk counterpart HH). The expression 
for T PPA is a generalization of the bulk counterpart |l(]] 
to the case of magnetic layers: a vanishing T pp is ob- 
tained for A = jlq in an agreement with the Mermin- 
Wagner theorem |l7] and small relativistic effects have 
to be considered in order to obtain a non- vanishing value 
of T PPA . The anisotropy energy A is taken here as an 
adjustable parameter. This is not a serious problem as 
the RPA Curie temperature has only a weak logarith- 
mic dependence upon A jlq ], and it is thus sufficient 
to know the order of magnitude of A. The latter is 
typically of the order of the dipolar energy 2irM 2 /V, 
where V is the atomic volume. In calculations we used 
A Co =0.052 mRy and A Fc =0.140 mRy. The evaluation 
of T PPA is facilitated by observation that it is propor- 
tional to the real part of the magnon Green function 
G m {z) — N^ 1 (z — i?(q||) _1 corresponding to a dis- 
persion law -E(qy) and evaluated at z — 0. The corre- 
sponding q|| -summation is performed for complex ener- 
gies from which the value at z = is obtained by an ana- 
lytic continuation technique jl8| . The sum for the evalua- 
tion of the spin- wave stiffness constant is non-convergent 
due to the RKKY character of magnetic interactions in 
metallic systems and to overcome this difficulty we have 
calculated it by a regularization procedure jl9| described 
in detail in Ref. 0. 

The calculated EEIs for the first 10 shells, magnetic 
moments, spin-wave stiffness constants, and the RPA and 
MFA Curie temperatures are summarized in Tables I and 
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FIG. 1. Magnon dispersion laws (left frame) and corre- 
sponding densities of states [in states/(meV- atom] for the 
Fe-layer embedded in fcc-Cu (full line), Fe-overlayer on 
fcc-Cu(OOl) (dashed line), and the free-standing Fe-layer 
(dashed-dotted line). We have set here A = 0. 

II for three limiting cases of magnetic Fe and Co layers, 
namely the free-standing fcc(OOl) layer, the overlayer on 
fcc-Cu(OOl) substrate, and the fcc(001) layer embedded 
in fcc-Cu host. 

Concerning the EEIs, the following general conclusions 
can be drawn: (i) A pronounced dependence of magnetic 
moments on the coordination number is found, namely 
their decrease with increasing number on nearest neigh- 
bors, the effect being stronger for the Fe layer; (ii) the 
EEIs are significantly enhanced (typically by a factor 2 
or more) as compared with their bulk counterparts; (iii) 
the EEIs depend strongly upon the presence of a sub- 
strate and a capping layer. The latter dependence is due 
to the RKKY character of the EEIs in metals: the cou- 
pling is not only mediated through the magnetic layer 
itself but also via the substrate and capping layer. This 
behavior is also clearly visible on the spin-wave spectra 
shown in Fig. 1 for Fe and on the exchange stiffness con- 
stant (Table II). We also present corresponding magnon 
densities of states (DOS) determined from the magnon 
Green function G m (z). A characteristic step of the height 
proportional to 1/D at the bottom of the magnon DOS 
accompanied by a pronounced van Hove singularity in 
the middle of the band are typical features of the two- 
dimensional bands with the nearest-neighbor interactions 
which are here only slightly modified by non-vanishing 
interactions in next shells. Interestingly, the spin-wave 
stiffness constants and Curie temperatures behave differ- 
ently as a function of the atomic coordination for Co and 
Fe layers, i.e., for cases of the free-standing layer, the 
overlayer, and the embedded layer. This behavior can 
be related to the values of leading EEIs in both cases, 
in particular to large antiferromagnetic couplings of 3rd 



FIG. 2. RPA Curie temperatures of Co overlayer on 
fcc-Cu(OOl) (full squares) and of fcc(001)-Co layer embedded 
in fcc-Cu (full circles) as a function of the magnetic anisotropy 
energy A. Note the logarithmic scale on the abscise. 

TABLE II. Calculated values of the spin- wave stiffness co- 
efficient D, T* PA , and T C MFA for Fe and Co magnetic layers: 
a free-standing (fs) layer, an overlayer (ov) on fcc(001)-Cu, 
and an embedded (em) layer in fcc-Cu. 
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and 4th nearest neighbors of Fe-based layers which effec- 
tively reduce the value of the spin-wave stiffness constant 
(see Eq. @), in particular for the free-standing layer. 
On the contrary, the Co-based EEIs have the prevail- 
ing ferromagnetic character giving thus increasing spin- 
wave stiffness constants due to increasing values of the 
EEIs with reduced atomic coordination. The antifer- 
romagnetic character of the EEIs for fcc-based Fe lay- 
ers is strongly enhanced as compared to the bcc-Fe case 
H while the prevailing ferromagnetic character of the 
EEIs for bulk fcc-Co || and for fcc-Co layers remains 
unchanged. 

The MFA Curie temperatures are typically of the same 
order magnitude as the corresponding bulk ones due to 
the fact that the reduced coordination is approximately 
compensated by the increase of the EEIs. This observa- 
tion is in a strong disagreement with experimental data: 
this failure is due to the fact that the MFA violates the 
Mermin- Wagner theorem due to the neglect of collective 
transverse fluctuations (spin-waves) and it is thus inap- 
propriate for two-dimensional systems. 

The RPA Curie temperatures as a function of the 
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FIG. 3. T c flpA of a Co (left) and Fe (right) overlayer on 
fcc-Cu(OOl) substrate covered by a cap layer of varying thick- 
ness. The dashed lines represent the embedded layer limit 
(infinite cap thickness) while the limit of zero cap thickness 
corresponds to the uncovered overlayer. 



the Co/Cu/Co trilayer with the varying thickness of the 
Cu-cap layer pl| . A similar behavior was also verified for 
the oscillatory dependences of T PPA and T c MFj4 which, in 
turn, are derived from the EEIs. It should be noted that 
amplitudes and phases of oscillations can be influenced 
by the thickness of magnetic layer and/or the presence 
of the disorder in the system. 

In conclusion, in view of the interpretation proposed 
here, the oscillatory behavior of the Curie temperature of 
Fe films as a function of the Cu-cap thickness as reported 
by Vollmer et al. jlj] would constitute the first direct ex- 
perimental evidence of the oscillatory RKKY character 
of exchange interactions in itinerant ferromagnets. 
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anisotropy energy A are shown in Fig. 2 for cases of Co 
overlayer on fcc-Cu (001) and fcc-Co(OOl) layer embed- 
ded in Cu. The weak logarithmic dependence of T PPA 
on A Q is obvious: T* PA varies by about 25% as A 
varies by an order of magnitude so that the results are 
not significantly influenced by our semi-empirical choice 
of A. The RPA Curie temperatures are strongly reduced 
as compared to the corresponding bulk values thereby im- 
proving on the MEA results. Nevertheless, they are still 
too large as compared to observed Curie temperatures 
of ferromagnetic monolayers (being of order 150 — 200 
K). It is unclear whether this is due to some inaccuracy 
of the theory or to some imperfections of the samples 
used in experiments. On the contrary, such important 
experimental facts as the strong influence of the metallic 
coverage on the Curie temperature [Q are well explained 
by our theory as illustrated in Fig. 3. The oscillatory 
character of T PPA around the value corresponding to an 
infinite cap, i.e., to the limit of the embedded layer, is 
clearly visible and it is in a qualitative agreement with 
the recent experiment of Vollmer et al. The origin 
of these oscillations can be traced back to the oscilla- 
tory behavior of the EEIs and it has the same origin as 
related oscillations of the interlayer exchange couplings 
found for the Co/Cu/Co(001) trilayer with a varying Cu 
cap-layer thickness [Q. These oscillations are due to 
quantum-well states in the Cu-cap layer formed between 
the vacuum and the magnetic layer which, in turn, influ- 
ence properties of the magnetic layer. We have verified 
that amplitudes of oscillations of the EEIs decay with 
the thickness d of the cap layer approximately as dr 2 . 
The same thickness dependence was also found for the 
related case of the interlayer exchange interactions for 
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